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ABSTRACT 


A study was made of Type I linear control systems, and cost functions, the integral 
of error squared, and the integral of error squared plus a weighting factor times the 
control effort squared. The effects of cost function minimization on characteristic root 
movement for specific third and fourth order systems are investigated to determine if 
they move into any recognizable patterns. The effects on the unit step response of the 
munimized systems are also determined. Performance of these minimized systems tn the 
presence of saturation is evaluated specifically to determine if the system response can 
be improved bevond that of function minimization alone. By minimizing the mtegral of 
error squared, IES, the characteristic roots do tend to the pattern of the IES standard 
form. Also, by weighting the control effort while minimizing the svstem error, the char- 
acteristic roots do tend toward the Butterworth pattern. Systems designed in this manner 


do perform better 1n the presence of saturation. 
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I. INTRODUCTION 


There are numerous methods which can be used to design linear control systems. 
The primary objective of the different techniques available is to give accurate control of 
the system output in steady state. However, the transient response of these systems re- 
presents the brunt of the problem for the design engineer. Of course the svstem must 
be stable, but more important it must reach steady state in an acceptable fashion. 
Hence, system performance criteria are given to the designer which dictate desired sys- 
tem performance during the transient and steady state portions of svstem operation. In 
general, these specifications, particularly maximum overshoot, settling time and steady 
state accuracy determine the ease or difficulty in designing a given system. In some cases 
they may make one method of design more suitable than another, but in all cases the 
dynamic response of a linear control svstem is determined by the roots of the charac- 
teristic equation. 

The classical methods for control design (1.e., Bode, root locus and parameter plane) 
give a dominant pair of complex roots and involve trial and error. More direct methods 
such as pole placement or the use of standard forms can significantly reduce the amount 
of time required to complete the design. But in all cases, no matter what design method 
is used, the roots of the characteristic equation are placed to give the desired system 
performance. Using the classical methods, the specific locations of the systems charac- 
teristic roots are not known. One simply designs a compensator, be it cascade or feed- 
back, and tests the system to see if it meets specifications. However, using the more 
direct methods, be it by standard form or pole placement, each of the characteristic roots 
are moved about specific locations in the s-plane until desired dynamic performance is 
achieved. 

For more complex svstems, it may not be reasonably possible to place all the char- 
acteristic roots in specific locations, or it may take an unrealistic number of trials to 
achieve the desired dynamic response using the trial and error classical methods. How 
then is this svstem to be synthesized ? One way is through the use of performance in- 
dices otherwise known as cost functions. Since the primary objective of most control 
systems is to give accurate control of the svstem output in steady state, if we minimize 
or control some function of the system error, E, the best possible svstem performance 


can be obtained. Some of the more common cost functions include 


J=| Evdt (1.0) 
Or 


J=| | Ejedt (1.1) 


Where the best system is defined to be that which provides the minimum yalue of the 
integral. 

If full state feedback is used. given that all states are available, each coefficient of 
the characteristic equation would be variable, and this would allow complete minimiza- 
tion of the given cost function. Hence. a particular closed loop root pattern would re- 
sult, producing a particular transient response which may or may not meet 
specifications. 

In order to gain a better appreciation of these cost functions a more common form 


given as 


Oo 


J= | (X'OX + U'RU)dt (1.2) 


Where X is defined as the state vector. The vector U is defined as the controlimom 
vector, and Q and R are symmetric matrices, called weighting matrices. Bv considering 
only the time variation of the states, the transient response, the cost function is reduced 


to 


J=| X'ONXd (1.3) 


Defining the states such that 





Me = lcs) 
eee Lal 
X,=-~.,... (1.6) 
: dt’ 


and giving zero weight to all derivatives, the cost function reduces to 


oO OO 


J=|) Xidt=| Edt (1.7) 


Better transient response can be obtained by weighting some or all of the other states. 
It must be pointed out, however, that the process of choosing values for the weighting 
matrix is arbitrary and the transient response obtained will depend on the values chosen 
for the weighting matrix. This in turn reduces the procedure to a trial and error method 
where the best design must be chosen bv inspection of the simulation results. [ Ref. 1] 
The main thrust of this study 1s to examine the relationships, if they exist, between 
specific cost functions and the closed loop roots of arbitrarily chosen plants. That is to 
say, for a given cost function, do the closed loop characteristic roots tend to a particular 
pattern? How does weighting the control effort of a given cost function affect the system 
root pattern? How does the system perform in the presence of saturation? Can anv 
specific conclusions or generalizations be reached? In order to answer these questions 


third and fourth order systems will be studied. 


Il. SYSTEM DEFINITIONS AND STANDARDIZATION 


The results obtained from the study of a particular control system are in most cases 
unique to that system. Performance requirements can make the design of one system 
easy to engineer in one case and difficult in another, even though the same performance 
requirements are used. In order to obtain data which can be compared, the methods 
used to derive this data will be common to all systems studied. These svstems, which 
will be defined later, are all pole plants and are arbitrarily chosen so that actual system 
performance is not known beforehand. 

This case study 1s concerned with the closed loop roots of subject systems. Even 
though the characteristic roots are not placed in specific locations as a matter of proce- 
dure, this is precisely what is done. A block diagram of a typical closed loop system is 
Shown in Fig. t. Here R is the system input, which will be a unit step function. E is the 
system error, U ts the control effort, and C is the system output. K is a scalar gain and 


G(s) 1s the open loop transfer function. 





Figure 1. Basic system diagram. 


Feedback control is ysed in this study as a matter of convenience. Simply put, 


feedback control is the use of a present system condition to influence its condition in the 


future. It is used here to shape the system’s transient response. Again it 1s pointed out 


that, in shaping the system's transient response, the characteristic roots are moved about 


the s-plane until a location 1s found which gives the desired response. The basic form 


of feedback control to be used is shown in Fig. 2. 


Figure 2. Feedback control arrangement: Unity feedback is preserved. 


Consider an all pole plant in the general form 


K 


G(s) a —| n—2 
aie aS +... sp Shame 


s*+A,_18 


R 


Using full state feedback, the feedback transfer function is 
Tis) sere ices ic, 
Then the characteristic equation 1s 
1 + G(s)H(s) =0 


so (A wt Kk 2, Js, oe take )s + (A, + Kk) = 0 


(2.3) 


(2.4) 


Since the roots of the, characteristic polynomial are a function of the polynomial’s 


coefficients, it is obvious that every root of the system can be located by adjusting the 


feedback coefficients. In doing this, the systems transient response is changed. There 


are many methods which can be used to adjust the feedback gains of these systems. 
Llowever, cost functions and function minimization will be used here. 

When one considers the use of cost functions in system design, this use can be 
thought of as a way of finding the best possible combination of system parameter values 
to minimize the value of the integral. It must be noted that the use of cost functions in 
system design does not guarantee that the design will meet specifications. Efowever, by 
minimizing the value of a given cost function, the best value for cach coefficient in the 
characteristic polynomial is specified. This in turn defines the locations or pattern for 
the characteristic roots. 

The characteristic root pattern for many design methods have been categorized into 
standard forms for the characteristic equation. Some of the more common standard 
forms have been used quite extensively in system design. They include the Binomial 
form, Butterworth form, and the Integral of time times absolute error form to name a 
few. Examples of these standard forms are shown in Fig. 3. These forms define the 
desired characteristic equations. State variable feedback can be used to obtain the re- 


quired coefficients. 





; Table fll. The Binomial Standard Forms 


eS ot i ees Ok eo 
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s+us 
57+ 2uns + wo! 
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Table IV. The Butterworth Standard Forms 





$+un 
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Table V. The Minimum ITAE Standard Forms, Zeto-Displacement-Error Systems 


= => = 
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Figure 3. Common characteristic equation standard forms: (From Ref. 1.) 


The cost functions considered in this study are 


OF 


o E’dt ( 


tJ 
Cay 
wae” 


and 


J=) (£°4+QUd (2.6) 
0 


where E is the system error, U is the control effort, Q 1s a weighting factor chosen by the 
designer based on the relative importance of E and U in the specific application. 
Algebraic evaluation of these integrals is a labor intensive procedure even for lower 
order systems. For higher order systems, the number of variables 1s further increased and 
algebraic evaluation of these integrals is a near impossibility. Evaluation of these and 
other cost functions can be completed readily by using computer simulation and func- 
tion minimuzation. This method of cost function evaluation 1s straight forward. The 


scheme is shown in Fig. 4. 







Variables 


FUNCTION MINIMIZATION Output 


R,E,U,C 


Figure 4. implementation of function minimization: (From Ref. |.) 


Dynanuc Siniulation Language, DSL, was used to implement this design method. 
This method works by using a program and a subroutine. The system to be simulated 
is inplemented in the program and variables are passed to the subroutine where the cost 
function is evaluated. This process is repeated iteratively until the cost function 1s min- 
imuzed. The operator must give initial values to start this procedure. While this method 
of system design is quite useful, one must be aware of constraints that its use places on 
the problem. It must be pointed out that the integral cannot be evaluated for an infinite 
period of time. Hence, a final time must be picked for the simulation. This does not 
present a real problem for the evaluation of the integral. One must simply understand 
that error will exist in the solution. The amount of error in the solution due to simu- 
lation time is a function of the plant under study and how much CPU time ts available. 
For most control systems, CPU time is not a problem. However, a compromise can 
usually be reached between CPU time and solution accuracy such that the error is in- 


significant. 


Additionally, in using this approach the parameters which are to be optimized must 
be independent and initial values for these parameters must be established. These initial 
Values must be picked such that the initial system is stable. If the initial svstem response 
Is not stable, the minimization routine mav not converge. Prudent choices for initial 
values will usually satisfy this criteria without difficulty. Other advantages of function 
minimization are 

1. Nonlinear cost functions can be evaluated. 


2. The nonlinear characteristics of the system are easily implemented in the simu- 
lation. 


3. Limitations such as amplifier saturation can be included in the simulation. 
4. Nonlinear compensators can be used and their parameters adjusted by the function 
minimization subroutine. 
By using function minimization a minimum, if desired, can be found for most cost 
functions. Even though some systems designed using this method will not be optimal, 


they will be the best svstem available given the allowable parameter adjustments. [ Ref. 


1] 


A. PLANT DEFINITIONS 

In an attempt to study systems which are some what complex, third and fourth or- 
der plants are considered. Second order plants are not considered because treatment of 
these svstems 1s quite common in design literature. Tvpe one systems are used because 
they give no steady state error to a unit step input. 

One of the problems which comes to light in a study of this nature is that for the 
systems which are to be studied, how can the data be generated so that 1t can be com- 
pared equitably? That is to savy, how can a comparison of apples and oranges be 
avoided? In order to establish a common Starting point, each of the plants are initially 
placed at the limit of stability. This forward gain value is decreased slightly if necessary 
in order to ensure that the function minimization subroutine will converge. 

1. PLANT 1 
The plant transfer function 1s 


K 


MOSSE IE + 6) oi 


Since we desire that the value of K be such that the plant is at the stability limit, the 


characteristic equation must be determined and the value of K found. The closed loop 


transfer function 1s 


G(s) 
Cel) <TH GG 


and since unity feedback is used 
H(s) = 1 
Therefore the cliaractemstic equatioms 


1 + G(s)H(s) =0 


K 


er Cea) 


= 0 


5) eee I25+ k =o 


(2.9) 


(2.10) 


(2.11) 


One 


To determine the required value for K substitute s =j@ into the characteristic equation 


to obtain 
— jo? — 8@" + lio + K=0 
Separating the imaginary parts of this equation and solving for w 
—w’ + 120 =0 
and 
@ = + 3.464 
Substituting this value of w into the remaining real part of equation (2.13) 
— 8(3.464) + K=0 


and 


The final plant transfer function is 


(2.13) 


(2.14) 


(2.15) 


(2.16) 


(2.17) 


96 


Cs s(s + 2)(s + 6) 


(2.18) 
The unit step response 1s shown in Fig. 5. This plant 1s clearly at the stability limit. The 
reason for testing the plant at this point serves two purposes. first, 1t verifies the cal- 


culations made so far, and second, it validates the computer simulation program. 


4] 
ra 


HUT PUT 


~ 
res 
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Figure 5. Third order system: Unit step response before cost function 1s mini- 
mized. Characteristic roots are located ats = +j 3.464. 


It 


Pee! oh (2 fel Oe 


The plant transfer function is 


K 


(OS TEENS + HEF) a 
Repeating the procedure used for plant 1, the characteristic equation is 
i (2.20) 
S(s + 2)(s + O)(s + 9) 
S lis qe sds-4 lOsceee=G 220 
Substituting s = j@ and separating the imaginary portion 
“Vig 10s (2.22) 
and 
w=+2.52 (2.23) 
Again substituting this value of w into the remaining real part of equation (2.21) 
(2.52) = Sap 52) ean (2.24) 
and 
hoe (225) 
The final plant transfer function is 
G(s) = we (2.26) 


The unit step response and closed loop root locations for this plant are shown in Fig.6. 
The reason for using this method to determine the value of K to place the svstem at the 
stability limit is two fold. First, it gives the gain value required, and second, it gives the 
value of wm where the jw axis is crossed. Again, it 1s obvious that this system is at the 
stability limit. The results obtained from simulating the system agree with calculations. 


Hence, as explained earlier, the calculations and the simulation program are validated. 


PET Pty] 


ua 


Figure 6. Fourth order system: Unit step response before plant is minimized. 


Characteristic roots are located ats = +) 2.52. 
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HI. THE INTEGRAL OF ERROR SQUARED 


If full state feedback is used to compensate the system, all coefficients of the char- 
acteristic equation become variable and a minimum of the integral of error squared ex- 
ists Onlv at @ =oo. To obtain a finite optimization the physical nature of the system 
requires that a constraint be placed on the bandwidth of the system. This is defined by 
specifving that the coefficient of the s° term in the characteristic equation should be w/ 
where \ 1S the Ondepoutie s.ctem: 

For this study, the specified systems are constrained so that w, is equal to 20 
rad/sec. This value was chosen as a matter of convenience. In order to do this one 
simply must keep in mind the relationship between the feedback gain k, and the plant 
forward gain value k, for the equivalent system. Mathematically, this relationship for 


the third order svstem 1s 
(Kok,) = (3.0) 


Because it is desired that the system frequency be 20 rad/sec, and the forward gain value 
kK, 1S equal to 96, one simply substitutes the known values and solves for the required 


value for k, 


96k) = $000 (3.1) 
ky = $3.333 (ao 


Repeating this procedure for the fourth order system we obtain 


(kok) = wo (3.3) 
493.1k, = 160000 (3.4) 
ky = 324.477 (3.5) 


This common value of @, is not chosen so that the third and fourth order systems can 
be compared. The two systems cannot be compared, and an attempt to do so would be 
like comparing apples and oranges. This value is used to help standardize each system 


while under study. 
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A. THIRD ORDER SYSTEM 

Carrving out the mininuzation procedure described earlier, it is found that the sys- 
tem’s unit step response is significantly iunproved. This response is shown in Fig. 7. In 
order to see how well cost function mininuzation works, the minimized response 1s 
compared to the unit step response obtained by implementing the standard form for the 
integral of error squared, IES. This response is also shown in Fig. 7. There is no dif- 
ference in the unit step response for the two design methods. The fact that there is no 
Steady state error 1s an expected result. In order to investigate the reasons for no dif- 
ference in the response for the two design methods, the characteristic equations and 


roots for both designs are shown in Table 1. 
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Figure 7. Third order system minimized response: Unit step response for mini- 


mized and IES standard form designs. 
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Table 1. Porson | HENS EQUATIONS AND ROOTS 
Chi wh AC LERAS TIC — | 
moO ALION Y2,3 
= 





s° + 20s* + 800s + 8000 =0 aiuaon 4.30 + j 26.14 


The results shown in Table | are expected since the systems unit step response are 


identical. Looking further into the results of the mininuzation process, an explanation 
is found that nught further explain the similarities in the unit step response. During the 
muninuzation process imitial guesses were made for the values of feedback gains required 
to minimize the cost function. After execution of the program. the values of the feed- 
back gains generated were iteratively placed in the minimization program. This process 
was repeated until the minimum was obtained. To see how well this worked, the final 
cost for the minimized and IES methods of design both equal 0.075. The simulation cost 
also equals 0.075. Each of the three values for cost 1s obtained by a different method. 
The minimization cost is the value determined from the function minimization program. 
The simulation cost is the value determined by placing the feedback gains generated in 
the function minimization program into a different simulation program. In this simu- 
lation program the cost function was again evaluated while generating the unit step re- 
sponse shown in Fig. 7. The integral of error squared cost value was also determined in 
the stmulation program, but it 1s the cost determined from evaluation of the cost func- 
tion for the standard form method of design. These values given for cost have no sig- 
nificance or real meaning. They simply represent the numerical value obtained from 
evaluation of the integral. Observing these cost, it is seen that the minimization and 
simulation cost are equal. The fact that they are equal serves to check the minimization 
and simulation programs for proper performance. The point is made by noting that the 
cost value for the IES method of design 1s equal to the cost for the minimization method 


of design. 


B. FOURTH ORDER SYSTEM 
The same basic steps carried out for the previous system are repeated for the fourth 


order system. It is understood that the actual results will not be the same, but the trends 


should be continued. As pointed. out earlier, the value of w, for this system is con- 
strained in order to control the system bandwidth. The value of &, which gives w, = 20 
rad sec was determined to be 324.477. Only the feedback gains &,, A, , and é, are allowed 
to be changed by the function mininuzation subroutine. Once the cost function ts min- 
imized and the final values for the feedback gains are determined, the minimized system 
is Simulated and the unit step response obtained. This minimized system is again com- 
pared to a system designed by using the standard form for the integral of error squared. 


The unit step response for these two designs are shown in Fig. 8. 
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Figure 8. Fourth order system minimized response: Unit step response for mini- 


mized and IES standard form design. 
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Again, the unit step response for these two methods of design are identical. The 
characteristic equation and closed loop roots for each method of design are also identi- 


Cal: 
s’ + 20s? + 1200s? + 1600s + 160000 = 0 (3.6) 


Roots are located ats = -2.10+)]31.l ands = -7.90 +) 10.14. On examination of the 


cost for these two designs it 1s found that they are again equal, 0.100. 
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IV. WEIGHTING THE CONTROL EFFORT 


In some practical applications, the number of adjustable parameters may be limited 
due to cost, weight, or space limitations. The design engineer must, taking these hmi- 
tations in to account, design the best possible system. As an example, if cascade com- 
pensation is used for the design, the compensator may be limited to one zero and one 
pole, where only the gain is adjustable. For higher order systems it may or may not be 


possible to determine an optimal svstem defined by a cost function such as 


f 


J= | (E°+QU)d1 (4.0) 


However, by using simulation and function minimization a minimum can be found for 
the cost function. It must be kept in mind that the system designed using this method 
will not be optimal, but will be the best available system given the parameters which are 
variable. This design may or may not meet specifications. [ Ref. 1] 

For the plants considered in this case study, how does the cost function affect the 
closed loop root locations? Does its use enhance or degrade system performance? 
Where do the closed loop roots of the characteristic equation tend as the weighting 
factor, Q, is changed? Specifically, Chang { Ref. 2] points out that as the weighting 
factor is decreased, the roots of the characteristic equation should tend to a Butterworth 
pattern. Shown in Fig. 9 is the Butterworth configuration for the characteristic roots. It 
is seen that the characteristic roots are equally spaced along the circumference of a circle 


whose radius is equal to @. 
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Figure 9. Butterworth pattern of characteristic roots: IHfere k is the system order. 
(From Ref 3.) 


A. PLANT } 


The open loop transfer function for this plant is given as 


96 


ON s(s + 2)(s + 6) 


(4.1) 
This system is implemented using the simulation and function minimization method 
introduced earlier. In this case the cost function is changed. Three basic steps are used 
to study this cost function. First, the system is simulated and the cost function is mini- 
mized with specific weights placed on the control effort, U. Second, the minimized sys- 
tems’ closed loop roots are determined and tabulated. Finally, the systems’ unit step 
response is obtained. This procedure is repeated with different weights placed on the 
control effort. 

Because the cost function will affect a given plant in a unique manner, there 1s no 


way of knowing if the system designed as a result of minimizing the cost function will 
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meet specifications. Further, because there is no method to determine how to choose 
values for the weighting factor, Q, it is normally chosen arbitrarily. HHowever, since it 1s 
known that, as the weighting factor is decreased, the characteristic roots should tend to 
a Butterworth pattern, the symmetric root locus is used to help make an educated guess 
for the proper value of the weighting factor, Q. 

For this svstem the symmetric root locus is deternuned by placing the characteristic 


equation in the form 


K? 


O° | 
GES) (s)( —s)(s + 2)\(s — 2)\(s + 6)(s — 6) Br) 





This symmetric root locus, shown in Fig. 10, was generated using the computer program 
EWALD. The symmetric root locus for this third order system shows how the root lo- 
cations radiate outward and tend asymptotically towards the Butterworth configuration. 
The process of minmuzing the cost function for various values of Q causes the roots of 
the characteristic equation to move outward along the locus to a specific value of @,. 
For each chosen value of Q, the minimum reached forces the characteristic roots into a 
closer approximation to the Butterworth pattern. It must be pornted out that as Q 1s 


decreased and the cost function is minimized, the value of w, is increased. 
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Figure 10. Symmetric root locus for third order system: Used to assist finding 


values for weighting factor, Q. 


Using this procedure, the true Butterworth pattern is reached only when @, = co. 
Shown in Fig.11 is the symmetric root locus with minimization results plotted. It is 
shown that the closed loop roots tend to a circular pattern. One of the unique features 
found during the minimization process was that for some values of the weighting factor, 
the characteristic roots did not plot on the symmetric root locus as expected. Efforts to 
determine why this happened were inconclusive. Flowever, this behavior may be attri- 


buted to reaching a local minimum during the minimization process. 
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Figure Il. Third order system minimization results: Characteristic root locations 
plotted on symmetric root locus shows root movement towards 


Butterworth pattern. 


Since the system characteristic roots are expected to tend towards a Butterworth 
pattern, the unit step response for the true Butterworth configuration is determined and 
compared to the minimization results. The true Butterworth unit step response was ob- 
tained by using the standard form for the third order system. This standard form was 
implemented using the same value of w, as those determined from minimizing the subject 


cost function. The Butterworth standard form for the third order system 1s 
s° + 2wos” + 2wgs + wh (4.3) 


The form of the system open toop transfer function ts 


ie 3 
Gs\= =t 


a (4.4) 
s+ 2wys” If 2ws 
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fie s+sterm Was simulated using w, equal to 15.96, 17.12, and 26.72 rad‘sec. sAs pointed 
out earlier, it is desired to make as true a comiparison as possible. Therefore, frequency, 
time of stmulation, and system type are the same. Figure 12 shows the results for each 
design when the weighting factor is equal to 0.0005. This establishes w, at 15.96 rad ‘sec. 
Comparing the minimized results to the standard form results it is seen that the systems’ 
unit step response are identical. Table 2 shows the characteristic equation and closed 
loop roots for Q = 0.0005 . 


MN 
(10 


f 

. 

_ 
f+ 


Figure 12. Third order system minimized unit step response: Comparison of 


Butterworth standard form design and minimized design, Q = 0.0005. 


US 


Table 2... THIRD ORDER SYSTEM DESIGN SUMMARY Q = 0.0005 


CHARA C Te Ras ike 


METHOD EQUATION 


BUTTERWORTH | 3+ 31.985? + 514.745 + 4124=0 | -15.92 
MINIMIZED — | s3 + 31.92s? + 509.575 + 4066.9 = 0| -15.96 | - 





As expected, it is shown that the characteristic equation and roots are essentially the 
same. Further investigation of these designs show that the cost incurred by minimizing 
the cost function is 0.104 for the Butterworth design and 0.125 for the minimized design. 
Only the integral of error squared was evaluated for the Butterworth design. The subject 
cost function was evaluated for the minimized response. This explains the difference in 
the two values obtained for cost. 

The second value of the weighting factor for which the subject cost function was 
minimized is Q = 0.0002. Shown in Fig. 13 1s the unit step response for each method 
of design. Again, it can be seen that the response for the two design methods are fairly 


close. Table 3 shows the characteristic equations and closed loop root locations. 
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Third order system minimized unit step response: 
= 0.0002. 


Figure 13. 
Butterworth standard form design and minimized design, Q 


THIRD ORDER SYSTEM DESIGN SUMMARY Q = 


Table 3. 
CEL ARACTERISTIC 
-17.12 -§.56 +] 14.83 


| BUTTERWORTH | s+ 34.245? + 586.195 + 5017.8 = 0 
ee L 16.19 | -7.36 +} 16.47 
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Figure 14. Third order system minimized unit step response: Comparison of 
Butterworth standard form design and minimized design, Q = 
0.000005. 
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Table 4. THIRD ORDER SYSTEM DESIGN SUMMARY Q = 0.000005 





Observing the characteristic equations and closed loop roots for this case, it is seen 
that significant differences exists. These differences are responsible for the differences in 
the unit step response for the two design methods for this case. The values of the cost 
functions are 0.062 for the Butterworth design and 0.060 for the minimization design 


method. 


B. PLANT 2 


Repeated here 1s the transfer function for the fourth order system. 


493.1 


CO —————E————— 4.5 
I) = 55 Dist OS 49) 2) 
The same basic steps carried out for the third order system are again completed for this 


system. The symmetric root locus is realized from 


tO 


+ rhe 


G(s)G( —s) = (4.6) 


(s)( —s)(s + 2)(s — 2)(s + O)(s — 6)(s + 9)(s — 9) 


Figure 15 shows the symmetric root locus for this system. 
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Figure 15. Symmetric root locus fourth order system: Used to assist in deter- 


muning values for weighting factor, Q. 


Again, the characteristic roots of this system move outward asymptotically towards a 
Butterworth configuration as the forward gain constant is increased by decreasing QV. 
As the cost function 1s minimized when Q is decreased, the Butterworth configuration 
is achieved as w, approaches infinity. For this case study w,) approaching infinity 1s the 
result of very light weighting on the control effort, U. One can obtain the Butterworth 
pattern for small values of w, by using pole placement. 

The results for the minimized system are shown in Fig. 16. As the control effort 
is decreased, the closed loop roots move outward along the root locus. Even though the 
closed loop root locations shown in Fig. 16 do not form a true Butterworth pattern, the 
roots do tend towards a very rough circular pattern. It must be pointed out that the 
process of minimizing the cost function is very tedious. Simply because a minimum 1s 
reached does not mere that the root locations will match the symmetric root locus 


when plotted. Numerous values of Q were cliosen and the system was mininuzed, but 
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all points did not match the root locus. An explanation for this behavior 1s not offered 


and is a subject for future study. 
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Figure 16. Fourth order system minimization results: Characteristic root to- 
cations plotted on symmetric root locus shows that roots do tend 


roughty to a circle. 


In order to gain further insight into the results obtained by minimizing the cost 
function, the unit step response for the true Butterworth configuration 1s determined and 
the two are compared. The true Butterworth unit step response was obtained by using 
the standard form for the fourth order system. This standard form was implemented 
using the same values of w, as those determined by minimizing the subject cost function. 


The Butterworth standard form for the fourth order system 1s 


5s + 2.6wys" os 3.40 + 2.6w9s + Wo : (4.7) 


b 
Since the subject plant is Type I, the system realized from the standard form is also Type 


[. The open loop plant transfer function is 


a2 


i 
= sys 775, oe) 
S + 2.6095 + 3.4daos° + 2.0m 95 
This system was simulated using w, equal to 11.15, 14.89, and 19.85 rad ‘sec. 

In order to complete as true a comparison as possible all system parameters which 
could be made equal are so. That 1s to say, frequency, time of simulation, and system 
tvpe are the same. The difference 1s the method in which these svstems are designed. 
Figure 17 shows the results for each design when the weighting factor, Q, is equal to 
0.001. This establishes an w, equal to 11.15 rad’sec. Comparing these two systems it 
is evident that the system which was minimized with the weighting factor Q = 0.001 1s 
a very close approximation to the Butterworth standard form. The transient portion of 
these two svstems 1s essentially over in approximately 0.75 seconds. The maximum 
overshoot of the Butterworth standard form system 1s slightly higher than that of the 
minimized system. FHlowever, taking these minor differences into account, the 
Butterworth standard form and the minimized systems’ step response are approximately 


ive same. 
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Figure 17. Fourth order system minimized unit step response: Unit step response 


for Butterworth standard form and minimized design. 


Earlier in this study, it was pointed out that the design of linear svstems is an exer- 
cise In positioning the roots of the characteristic equation. Even though this is not done 
here as a matter of procedure, this is precisely what is achieved by evaluating the cost 
function or using the standard form to realize a particular design. Where do the char- 


acteristic roots go for these systems? Table 5 shows the results for these systems. 
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Table 3. CHARACTERISTIC ROOT LOCATIONS Q = 0.001 


METHOD hi 


BUTTERWORTH “234 1031 
Piles VIEZED 4.) fae | 913 -11.56 + 3 4.08 





From Table 5 the root locations for the Butterworth design and the minimized design 
are nearly the same. This further explains the similarity in the unit step response for 
these two designs. 

The second value of the weighting factor for which the subject cost function was 
minimized is Q = 0.0001. At this weighting, the frequency of the system was found to 
be 14.89 rad/sec. Again using the Butterworth standard form, the characteristic 


equation for this system 1s 
s* + 38.7145? + 753.8215? + 8583.362s + 49156.255 =0 (4.9) 


The characteristic equation which was determined from the minimization of the subject 


cost function is 
5° + 41.355° + 799.53s° + 8844.96s + 49054.96 =0 (4.10) 


Note the values of the coefficients for the Butterworth characteristic equation, equation 
(4.9) and the minimized svstem characteristic equation, equation (4.10). The unit step 
response for these two designs are shown in Fig. 18. The same characteristics stated 
earlier for the previous design are again true with one exception. The systems overall 
response is faster, as expected, because of the increased value of w,. Table 6 shows the 


characteristic root locations for these systems. 
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Figure 18. 
systems unit step response Weighting factor Q 


CHARACTERISTIC ROOT LOCATIONS Q = 0.0001 


Table 6. 
[wera g e 
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Again, the characteristic root locations for the Butterworth and the minimized designs 
@resvery close. I his explains the similarity in the unit step response for these two svs- 
foms. 

The unit step response for the remaining weighting factor Q = 0.00001 1s shown in 
Fig. 19. The trends established earlier are maintained and agree with expected results. 


Table 7 shows the characteristic root locations for these svstems. 
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Figure 19. Fourth order minimized system: Butterworth, and minimized svstems 


unit step response. Weighting factor Q = 0.00001. 
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Table 7. CHARACTERISTIC ROOT LOCATIONS Q = 0.00001 


METHOD see 
BUTTERWORTH 7 Sone nies hehe ol 
MINIMIZED Wis ye aie ty oe 






Viewing the results shown in Table 7, 1t 1s obvious why the unit step response for these 
two design methods are in such close agreement. For all practical purposes, these results 


can be considered equal. 
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V. SYSTEM PERFORMANCE WITH SATURATION 


For most practical linear control systems, there are limits on the linear operating 
range. These limits are usually given in the desired performance specifications for each 
system. In some cases, these limitations place severe constraints on the system to be 
designed, especially the dollar cost and physical weight of the desired system. [n most 
cases, the size of the power source for the desired system will determine how well the 
system performs when pushed to extremes. How well a system performs when in satu- 
ration is of particular concern. Large inputs can be damaging to mechanical compo- 
nents for control svstems. These inputs can cause electrical components to fail. For 
these reasons most control systems are designed such that they will saturate at lower 
input power levels. Friedland [ Ref. 5] points out that the effect of occasional control 
saturation 1s usually not verv serious, and systems which never saturate are probably 
over-designed. However, if the system almost always saturates it probably will not per- 
form satisfactorily. Of course, bang-bang control svstems are not included here. The 


cost function considered for saturation study 1s 


io) 


J=) (E°+QU%)d (5.0) 


A. MINIMIZED SYSTEM SATURATION RESPONSE 

In order to determine the system performance in the presence of saturation a limiter 
is placed in the system as shown in Fig. 20. For this system, it is desired to limit the 
control effort. The procedure used to determine how well the system performs when in 
saturation is to note the unit step response for the system when the plant is not in sat- 
uration, just into saturation, and well into saturation. The limiter levels which give these 
results are not particularly important. It 1s desired to obtain the overall system per- 


formance. 
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Figure 20. Saturation study system block diagram: Simulation method used to 


study system saturation performance. 


1. Third order system 

Figure 21 shows the unit step response for the minimized system when in sat- 
uration. Output MIN represents the system response without saturation. Output 
MINA represents the systems response with a small amount of saturation, and output 
MINB represents the system response when Well into saturation. It 1s known that this 
or most other systems will not perform satisfactorily when allowed to saturate heavily. 
This system was placed far enough into saturation to illustrate how the system's per- 
formance begins to deteriorate. Clearly the system performs satisfactorily with a small 
amount of saturation, even though the response is slower. However, when the system 
is driven into heavy saturation the system’s performance begins to deteriorate. Note the 
marked increase in the peak overshoot for this case. The system is beginning to go so 
far into saturation that it cannot respond to changes in the control input. To illustrate 


this point, Fig. 22 shows the output of the limiter. 
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Figure 21. Third order system saturation step response: Outputs MIN, MINA, 
and MINB represent the system response with no saturation, mild sat- 


uration, and heavy saturation respectively. 


Here U 1s the svstem control effort previously defined. LIMUA and LIMUB are limiter 
outputs for mild and heavy saturation respectively. These results indicate for mild sat- 
uration that the control effort 1s limited on the peak of the positive going portion of the 
control effort only, as seen in output LIMUA. When the system is driven into heavy 
saturation the control effort is limited on both the positive and the negative-going peaks. 
In each case where the svstem 1s allowed to saturate, the control effort to the system is 
less than the normal no-saturation value. Thus. when in saturation, the svstem receives 


power at a lower rate and cannot respond as fast as the unsaturated system. 
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Figure 22. Third order system control effort changes: Variation of control effort, 
U and hmiter output with mild saturation, LIMUA, and heavy satu- 
ration, LIMUB. 


If minimizing the cost function, equation (5.0), 1s viewed as determining the area 
under the error and control effort curves, when the svstem 1s allowed to saturate this 
area increases. Ilence, the cost of operating the system should increase. This point 1s 
verified by evaluating the cost function as the saturation level 1s increased. Shown in 


Table 8 are the results obtained. 
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Table 8. VALUE OF COST AT DIFFERENT SATURATION LEVELS 











The results seen in Table 8 illustrate that the value of the cost function does increase as 
the saturation levels are increased. Hence, the area beneath the error and control effort 
curves 1S increased. 

Chang [ Ref. 2 ] proposed for the subject cost function, that the system per- 
formance can be improved when the system is allowed to saturate at low levels. This 
point is verified by viewing the system’s step response shown in Fig. 21. Observe in the 
system output for the mild saturation case, MINA, the peak overshoot is slightly less 
than that for the unsaturated case, MIN. Additionally, the speed of svstem’s response 
iemrecuced, However, this is expected. Whether a slower system with less overshoot is 
more desirable than a faster system with more overshoot is a question to be answered 
by considering the design specifications. 

2. Fourth Order System 

The same procedure used for the third order system is repeated for the fourth 
order system. Also, the same cost function, equation (5.0), is considered. Shown if Fig. 
23 1s the unit step response for this system. Outputs MIN, MINA, MINB represent the 
system performance for no saturation, mild saturation, and heavy saturation respec- 
tively. Even though the third and fourth order systems performance cannot be com- 
pared, the trends found for the third order system are also displaved by the fourth order 
system. Viewing Fig. 23 for the mild saturation case, output MINA, the system’s re- 
sponse is slower and has a lower peak overshoot than the no saturation case, output 
MIN. For the heavy saturation case the system response is even slower than the mild 
saturation case. However, note that the peak overshoot 1s higher. If the svstem is driven 
further into saturation than that shown for the heavv saturation case, the svstem output 


becomes oscillatory. 
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Figure Fourth order system saturation response: Outputs MIN, MINA, and 
MINB represent the systems’ response with no saturation, mild satu- 


ration, and heavy saturation respectively. 


To further illustrate how the system performance is affected by saturation, the 
limiter output is shown in Fig. 24. System outputs U, LIMUA, and LIMUB represent 
the system’s performance for no saturation, mild saturation. and heavy saturation re- 
spectively. Again it is shown how saturation affects the system control input. For the 
mild saturation case the control input is limited only during the highest positive peak. 
For the heavy saturation case, the control input is limited on both positive and negative 
peaks. Clearly, the svstem’s contro! effort is limited to illustrate that the svstem 1s op- 


crating at maximum capacity. Hence the reason for the system outputs shown in Fig. 
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23 1s illustrated. When the svstem is in saturation, output is produced at a slower rate 
than when no saturation 1s present. 
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Figure 24. Fourth order system control effort changes: Variation of control ef- 


fort. LU and limiter output with mild saturation. LIM(UA, and heavy 
saturation, LIMUB. 


By evaluating the value of the cost function for the saturation levels shown in 
Fig. 24, it is expected that the cost will increase as the system is driven further into sat- 
uration. These results are shown in Table 9. These results show that the value of the 


cost function does indeed increase as the system 1s driven further into saturation. 
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Table 9. WALUE OF COST AT DIFFERENT SATURATION LEVELS 











In order to determine if the system performance 1s improved when operated with 
mild saturation, the design specifications must be considered. Figure 23 shows that for 
the case of mild saturation the system does have a lower peak overshoot but its response 
time 1s increased. A determination of whether or not this response is more desirable than 
the no saturation case, faster but with more overshoot, is based on the desired system 


performance specifications. 


B. COMPARISON WITH BODE DESIGN 

In order to further investigate svstem performance in the presence of saturation, the 
system Which was designed by weighting the control effort is compared to a Bode design. 
This comparison 1s completed for one specific reason. Chang [ Ref. 2] states that a 
svstem design by weighting the control effort is less affected by saturation. This effect 
is investigated and verified. 

Design of the Bode system is completed by simply designing a cascade compensator, 
or filter to give the desired output. In this case, the desired output is that which 1s sim- 
ilar to the output of the weighted minimized design. Only two specific designs are 1m- 
plemented. one each for the third and fourth order systems. For each of these systems 
only one value of the weighting factor was chosen to be designed. The third order Bode 
system was designed to approximate the minimized design for Q = 0.005, and the fourth 
order system was designed to approximate the minimized design for Q = 0.001. A block 


diagram of this implementation 1s shown in Fig. 25. 
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Filter Plant 


Figure 25. Implementation of Bode design: System arrangement for Bode design 


study. Note that the control effort, U, is generated by the filter. 


1. Third order system 


The transfer function for the filter of this design is 


pe OSG 2a ©) - 
We) "(5 + 20)(5 + 50) ) 


The plant transfer function ts 
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Upon completion of the design, the system was stmulated and the output compared to 
that of the weighted minimized design. Shown in Fig. 26 is the unit step response for 
the Bode and minimized designs. Even though the response of the two designs are not 
identical they are close enough to illustrate the main objective. Figure 27 shows the 


control effort variations for both design methods. 
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Figure 26. | Bode and minimized design unit step response: Outputs BODE and 


MIN represent Bode and minimized response respectively. 


While the overall shape of the control effort for both designs are basically the 
same, there are significant differences. Note that for the minimized response, UM, the 
control efforts starts at a maximum value and goes to zero. For the Bode design how- 
ever, the control effort, UB, starts at zero, increases to a maximum. and decay¥s to zeman 


These differences are expected because of the two different designs used. 
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Figure 27. Control effort for Bode and minimized design,no saturation: UB and 
UM represent control effort for the Bode and minimized designs re- 


SWECHIVEl.. 


Since it is desired to determine how saturation affects the weighted minimized 
design, the Bode design is used as an alternative design for comparison purposes. In 
order to introduce saturation into the Bode design, Figure 28 presents a block diagram 
illustrating limiter placement. Limiter placement for the minimized system 1s shown in 
fice 20. Three different levels of saturation were introduced: + 10 V, + 5 V, and + 2 


V. A unit step input was applied to both systems for all saturation cases studied. 
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Figure 28. —_ Bode design limiter placement: Limiter placed at plant input to sim- 


ulate saturation. 


Figure 29 presents the unit step response for these two designs with saturation 
set ata +10 V level. At this saturation level there is no noticeable change in outputs for 
either designs. Figure 30 shows the changes in the control effort. Note that the Bode 
design control effort, LBUA, is not yet affected by saturation. Therefore, its unit step 
response is unchanged. Ilowever, saturation does effect the minimized system's control 


effort, LMUA. However, there is no noticeable change in the system's step response. 
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reure 29. Bode and minimized designs, saturation level + 10 V: Output 


BODEA, and MINA represent the Bode and minimized designs re- 


SPeEcuvely: 


Fu 


ful 
1 
on 7 — 
ee A Te 
ee eS 
Ta 


1 MILA 


T UETTA 


tt 
ia 


Figure 30. Control effort changes, saturation level + 10 V: Outputs LBUA, and 


LMUA represent lnuter output for the Bode and minimized designs. 


Placing the system further into saturation is accomplished by reducing the 
limiter levels. Shown in Fig. 3] ts the unit step response for the systems with saturation 
set at + 5 V. Clearly, there are noticable differences in the output for the minimuged 
design. However, the changes in the Bode design are more drastic. Again, the svstem 
outputs are BODEA and MINA for the Bode and the nunimized designs respectively. 
Note that. while the mininuzed design has shghtly more overshoot, its overall shape 1s 
maintained. The Bode design, however, shows significant changes. Clearlv the over- 
shoot is no longer present and the svstem has become critically damped. Variations in 


the Systems control ellont are snowm Wee. 


BIA 


MM | ty 


Figure 31. Bode and minimized system response, saturation level + 5 V: Outputs 
BODEA and MINA represent Bode and minimized design response 


respectively. 


System outputs are as previously noted. Clearly, the effects of the limiter is illustrated. 
For the minimized system saturation takes place on positive and negative peaks of the 
control effort. Saturation takes place on only the positive peaks however for the Bode 
design. Note the difference in the response as seen in Fig. 31. The Bode design is more 


significantly degraded than the minimized design for the same saturation. 
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Figure 32. Control effort changes, saturation level + 5 V: Outputs LBUA and 


LMUA represent Bode and minimized designs respectively. 


In order to see if this trend is continued, Fig. 33 shows the unit step response 
for these systems when saturation is set at +2 V. Clearlv, there are significant differences 
in the response for both designs. The Bode design 1s clearly over damped. The mini- 
mized response has maintained its original shape even though its peak overshoot 1s sig- 
nificantly higher. Figure 34 shows the changes in the control effort for these two 


designs. It is obvious that the minimized system 1s further into saturation than the Bode 
design. 
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Figure 33. Bode and minimized system response, saturation level + 2 V: Outputs 
BODEA and MINA represent Bode and minimized designs respec- 


ele 


Note output LMUA; it swings from positive to negative and back into positive satu- 
ration. The Bode design however, output LBUA, 1s only saturated on the positive peak. 
Based on the amount of time each svstem 1s In saturation, the minimized svstem re- 
sponse is better than the Bode design. In order to gain further appreciation of how sat- 
uration effects these two designs, Table 10 presents the cost for each design at each 


Saturation level. 
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Figure 34. Control effort changes, saturation level + 2 V: Output LBUA and 


LMUA represent Bode and minimized designs respectively. 
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Table 10. SATURATION COST SUMMARY BODE AND MINIMIZED DE- 
SIGNS 


SATURATION LEVEL | BODE DESIGN COST cost DESIGN 


= 10 
: 






Here the weighted cost function shown in equation (5.0) is evaluated. Viewing 
these results it is clear that as each system is driven farther into saturation the cost in- 
creases. Even though there is a more significant cost increase for the minimized design, 
its response 1s clearly more acceptable than the Bode design. If the cost 1s computed 
onlv for the integral of error squared, IES, for these two designs totally discounting the 
effects of weighting the control effort, these trends are continued. Table 11 shows these 
Sectcm Note that for each level of saturation the cost mcreases, but the cost for the 
minimized system is less than that for the Bode system in all cases. Hence, the effects 


Olesatiration are less. 


Table 11. IES SATURATION COST SUMMARY BODE AND MINIMIZED 
DESIGNS 


VISE IZED DESIGN 


eePURATION LEVEL BODE DiEslGN, COsr COST 





2. Fourth order system 


Repeating this process for the fourth-order system, the trends established earlier 
are clearly evident and maintained. The transfer function for the final filter design for 


this svstem 1S 


{sy 
~} 


_ 811.64(s + 2)(5 + 6)(s + 9) 


SS) (s + 15)(s + 20)(s + 200) Se) 


The plant transfer function 1s 


493.1 


CS) — ean EO 


(5.4) 
Again the system is simulated and the response of the Bode design is compared to that 
of the weighted minimized design response. Figure 35 presents the unit step response 
for each design method. Clearly the response for these two designs are close enough for 
the saturation study to be conducted. The Bode design has more overshoot, but this will 
help illustrate the desired effects. Also shown in Fig. 36 is the changes in the control ef- 


fort for these two designs. 
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Figure 35. Bode and minimized design unit step response: Output for these de- 


signs are reasonably close; no saturation 1s present. 


Clearly the shapes of the transient portion of the response are the same with one ex- 
ception. The control effort for the Bode design, output UB, starts at zero builds to a 
maximum and proceeds to zero. However, the minimized design starts at a maximum 
value and proceeds to zero. These differences are due to the poles and zeros of the filter 
transfer function for the Bode design. Even though these differences exist in the control 


effort for the two designs, the output for each 1s as desired. 
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Figure 36. Control effort for Bode and minimized designs, no saturation: Outputs 


UB and UM represent Bode and weighted minimized designs respec- 


tively. 


Shown in Fig. 37 is the unit step response for these designs when saturation 1s 
introduced at a level of +10 V. Clearly the effects of saturation are obvious. iNotewmie 
drastic change in the Bode response. Even though it finally reaches steady state, it 1s 


clearly over damped. Figure 38 presents the control effort for these two designs. 
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Figure 37. Bode and minimized design step response, saturation level + 10 V: 
Effects of saturation are obvious. Note how significantly the Bode re- 


sponse, output BODEA, is affected. The minimized response remains 


unchanged. 


By limiting the maximum excursion of the control effort to 10 V, it is never allowed to 


reach its maximum Value. Saturation clearly has a significant effect on the Bode design. 
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Figure 38. Control effort changes saturation level + 10 V: Saturation at this 


level does not allow the control effort for the Bode design, output 


LBUAtosreach its Wiad aie: 


Bv allowing this system further into saturation the Bode response 1s further 
deteriorated. Shown in Fig. 39 is the step response for a saturation level of 4; 5V mac 
expected the Bode response 1s further degraded. Note that the minimized design response 
maintains its shape with only a minor increase in maximum overshoot. The control ef- 
fort variations for this level of saturation are shown in Fig. 40. 
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Figure 39. Bode and minimized design step response, saturation level + 5 V: 
Bode response, output BODEA, significantly affected. Minimized re- 
sponse maintains its shape with only a minor increase in maximum 


Overshnoot 


Note that the control effort for the minimized design, output LMUA, is in saturation 
on both positive and negative peaks. The effect on the step response 1s minimal. How- 
ever, the control effort for the Bode design 1s saturated onlv on its positive peak, output 


LBUA, and for less time, but its effect on the step response 1s clearly more significant. 
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Figure 40. Control effort changes, saturation level + 5 V: Control effort for 


Bode design, LBUA, is in saturation for onlv a short period, but its ef- 
fect 1s significant. 


The unit step response for a saturation level of + 2 V is shown in Fig. 41. Again 
the Bode design is significantly affected by saturation. But note the changes in the min- 
imized design response, output MINA. Even though there 1s more overshoot, this design 


response 1s clearly more acceptable. Figure 42 shows the changes in control effort for + 
2 V level of saturation. 
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Figure 41. Bode and minimized design step response, saturation level + 2 V: 


Bode design 1s severely affected by saturation, output BODEA. Miuni- 


nuzed design response deteriorated but less than Bode design. 


Note that for the first time the control effort for the Bode design is saturated on both 
positive and negative peaks, output LBUA. Also note that the minimized design is lim- 
ited on positive peaks for considerably more time, output LMUA, but the effect is not 


as significant as for the Bode design. 
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Figure 42. Control effort changes, saturation level + 2 V: Bode design signif- 
icantly affected by saturation, output LBUA. Minimized design 1s 


clearly deteriorated but not as significantly as Bode the design. 


Clearly the minimized design response is more acceptable even though its response 1s 
slower than the no saturation case. The speed of the response is slower, as expected, 
because each svstem receives power at a lower rate. 

In order to determine how saturation affects the cost for each design while in 


saturation, Table 12 presents the cost for each design at each saturation level. 
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Table 12. SATURATION COST SUMMARY BODE AND MINIMIZED DE- 
SIGNS 


SMIERATIONTEVEL | BODEDESIGN cost | ZED DESIGN 


0.246 
0.243 
0.29 
0.397 








Note that for each level of saturation there is a cost increase. However, the cost for the 


minimized system is considerably less than the cost for the Bode designed system. If cost 
is evaluated for only the integral of error squared, IES, for both design methods, totallv 
discounting the effect of weighting the control effort for the minimized design, these 


trends are repeated. Table 13 shows these results. 


Table 13. IES SATURATION COST SUMMARY BODE AND MINIMIZED 
DESIGNS 


SATURATION LEVEL | BODE DESIGN COST — cost 


NONE 










As expected, these results show that as the amount of saturation 1s increased, the cost 


is also increased. By evaluating the IES cost for both design methods, the difference in 
cost for the two design methods are more pronounced. However, only the minimized 
design was specifically designed for minimum cost. The Bode design was completed with 
no consideration given to cost. 

Based on these results for both the third and fourth order svstems, the mini- 
nuzed design response is clearly superior to that of the Bode design when operated in the 


presence of saturation. 


Vir CONCEUSIONS 


As stated in the introduction, the purpose of this study was to examine the re- 


lationships, if they exist, between the cost functions 


oo 


J=| Evdt (6.0) 


and 


oo 


J=| (22+ Q0U%)ar (6.1) 


and the closed loop roots of arbitrarily chosen Type [| all-pole plants. First, a determi- 


nation was made on how the cost functions affected the characteristic roots for the 


plants 
96 
(OSTEND) “J 
and 
' 493.1 
oO s(s + 2)(5 + 6)(s + 9) ce 


Specifically, it was desired to determine if the characteristic roots tend to a particular 
pattern as a result of minimizing these cost functions. Second, the effect of weighting 
the control effort on the characteristic roots was considered in order to determine if the 
characteristic roots again formed a unique pattern. Finally, the effect on system per- 
formance in the presence of saturation was considered to determine if the systems’ re- 
sponse could be improved. 

As a result of minimizing the cost function shown in equation (6.0), the integral of 
the error squared, the characteristic roots do tend to the same locations found by using 


the standard form for the integral of error squared. However, in order to reach this re- 


68 


sult, using function minimization, a specific technique must be used. First, the desired 
system bandwidth must be known. This value will usually be dictated by the given sys- 
tem specifications. Second, special consideration must be given to the initial values used 
in the function minimization program. If they are chosen arbitrarily, the minimization 
program may converge to a local minimum, or it mav not converge at all. After nu- 
merous attempts were made at choosing initial values using many different methods, all 
of which failed to give satisfactory results, a successful method was finally determined. 
Since no prior knowledge of system behavior was known, the standard form for the in- 
tegral of error squared was used to help pick these initial values successfully. 

To illustrate this procedure, the standard form for the third order integral of error 


squared system Is 
Sn ea ea. (6.4) 


For this study the value for @, was chosen arbitrarily. This value, 20 rad‘sec, was sub- 


stituted into the standard form giving 
s° + 20s* + 800s + 8000 (6.5) 


Since full state feedback was used to implement the function minimization program, the 


characteristic equation for this system was placed in the form 
s° + (S + 96k)s* + (12 + 96K,)s + 96K (6.6) 


Equating coefficients, the feedback gains are 


ky = 83.33 (6.7) 
k, = 8.21 (6.8) 
ky = 0.125 (6.9) 


Using these values as initial guesses, the function minimization program converges rap- 
idlv. Finally, the amount of time chosen to run the simulation must be carefully deter- 
mined. If it is not, considerable time can be spent minimizing the cost function and, 
hence, the dollar cost for computer CPU time is increased significantly. System per- 
formance specifications should help determine a good choice for simulation run time. 
As a consequence of these results, for linear control systems in which it 1s desired 


to minimize system error. the standard form can be used to complete the design. If this 
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design meets specifications no further work 1s required. However, if this design does not 
meet specifications other design methods must be used. 

As a result of weighting the control effort, the characteristic roots were found to 
tend to the Butterworth pattern as pointed out by Chang [ Ref. 2 ]. However, there is 
no technique for choosing the weighting factor, Q. The symmetric root locus was used 
in this study to help determine values for the weighting factor. However, its use was 
limited and only served to indicate that the system characteristic roots were moving in 
the proper direction. For some choices of the weighting factor the characteristic roots 
did not always agree with the root locus as expected. The actual reason for this behavior 
is unknown. This behavior may indicate that a local minimum was determined. This 
behavior 1s recommended for further study. 

As a result of evaluating the system performance in the presence of saturation it was 
determined that the system performance may be improved. If the system is allowed to 
saturate moderately, the performance is not significantly affected. Even though the 
system response time is increased, it 1s usually acceptable. However, if the system is 
heavily saturated, the performance rapidly deteriorates and is unsatisfactory. 

The comparison of Bode and minimized designs gave interesting results. These re- 
sults show that the system designed by weighting the control effort when minimizing the 
cost function does perform better in the presence of saturation. However, it must be 
noted that the Bode design was completed using cascade compensation, and the mini- 
mized design was completed using state variable feedback. The fact that different 
methods were used to implement each design may have affected the results. However, 
the results obtained do support Chang’s proposal [ Ref. 2] that systems designed by 
mininuzing the cost function in which the control effort 1s weighted, equation (6.1), 
perform better in the presence of saturation. Comparison of the minimized design with 
svstems designed using other methods is recommended for future study. 

Saturation studies for the integral of error squared were not conducted and are re- 
commended for future study. Also, the saturation response for the systems studied was 
limited to a single value of the weighting factor in order to keep the study manageable. 
It is understood that generalizations should not be made based on insufficient informa- 
tion. It is not put forth that these conclusions apply to all hnear systems. These results 
do indicate trends in system performance. Before they can be applied categorically to all 


linear control systems, they do require and are recommended for further study. 
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APPENDIA 


The following programs were used to produce the data used in this case study. These 
programs were written for use with Dynamic Simulation Language/VS, DSL. Other 
programs used but not listed were ALCON and EWALD/VS. 


a alawleulsalsule seftuwanton' sams aefae'a jan! ewan ewen se owes eee ewantow'a owawe ewawe oe) te Cheba 
[CCA GR ise le hoe be ke Ae bie Gi i oie ie Gi Oi oie iS oe i ii i i i a i ee ek Re ek ok ete ie ok ei 


% THESIS RESEARCH : FUNCTION MINIMIZATION 3RD ORD SYS '- * 
* THIS PROGRAM MINIMIZES THE COST FUNCTION, THE INTEGRAL* 
* OF ERROR SQUARED AND RETURNS FEEDBACK GAINS REQUIRED. * 


ve fies VARIABLE KO IS CONSTRAINED IN THIS PROBLEM w 
“ eo) 1HAL WO IS FIXED. ¥ 
Su i Se SERS Arr Eig hear rt rte tiie hi be cirided arbiter citirbir titer artis third sei ke tortie tee ark te rected seer ied arts 
D COMMON/HANDJ/FLAG, COST ,K1,K2 


ieee MINIMIZATION OF INTGRLCE**2),3RD ORD SYSTEM KO FIXED 


feet?) = Kl, X(3) = K2 

INCON ICO=0. 0 

Semol Ki0=1. 9583333, K20=. 02083333 
CONST KO = 10. 41677 


INITIAL 
WAGEICAG. GE. 0.0)G0O TO 5 
Kl = K10 
K2 = K20 


5) CONTINUE 
FLAG = FLAG + 1 


ate 
ray 


ale 
¢s 


DERIVATIVE 
no) Oe oS TEE 0. 0) 
E> Re=~OUTPDT 
i= hOee Re Klee X20 
OR ea 
Move 96.2 = 6 KID = 12% 5D 


B20 = INTGREC ICO; X3D) 
XD = INTGRLC ICO, X2D) 
OUTPUT = INTGRLCICO,xXD) 
FCNER = E**2 
El = INTGRLCICO,FCNER) 
CONTROL FINTIM = 5.0 
TERMINAL 
COST = El 
END 
STOP 
wv 


FORTRAN 
“= INITIALIZATION OF PARAMETERS FOR HOOKE SUBROUTINE 


IMPLICIT REAL*8(A-H,0-Z) 

DIMENSION X(2),STEP(2),Q(2) ,QQ(2) ,W(2) 
STEP(1) =0. 001 

STEP(2) = .0001 


N=2 
ITMAX = 500 
CFTOL = .0001 
ALPHA = 2.0 
BETA = 0.5 
IPRINT = 0 
MINMAX = -1 
CALL HOOKE(X,STEP,N,ITMAX,CFTOL, ALPHA, BETA,CF,Q,Q0,W, 
LIPRINT , MINMAX) 
STOP 
END 


Tevcredevevedevedededvevededede de tedeve sede desede seve dete seve te Tete ne dove ke seve dete te desevede te deve He fede de ved Fe ve ve de 
% THESIS RESEARCH : FUNCTION MINIMIZATION 4TH ORD SYS ¥ 
* THIS PROGRAM MINIMIZES THE COST FUNCTION, THE INTEGRAL* 
* OF ERROR SQUARED AND RETURNS FEEDBACK GAINS REQUIRED. * 
¥ KO IS CONSTRAINED TO LIMIT WO. 7 
Poe ese ab ok Ok i ake ak Ok ak i oi or Oo dele el kif Se Re OS SR ee ie) ley 
Phy 

D COMMON/HANDJ/FLAG , COST ,K1,K2,K3 

TITLE MINIMIZATION OF INTGRLCE**2),4TH ORD SYSTEM, KO FIXED 
Phy 


fee) — Ki, X(3) = K2;X(4) = K3 
INCON ICO=0. ,K10=32. 2287568,K20=2. 2632000, K30=6. 0839586E-03 
CONST KO = 324.477793 


INITIAL 
PCr uAG: Gk. 0.0)GO TO 5 
Kl = K10 
K2 = K20 
K3 = K30 


5 CONTINUE 
FLAG = FLAG + 1 
* 
se 


DERIVATIVE 
he IO eo ln, (0,0) 
be fh — OUTRUT 
M1 = KO *“E - K3 * X3D 
I = M1-K2*X2D 
Ce= 1 = Ke aD 
¥ 
X4D =493. 1*U-17*X3D-84*X2D-108*XD 
X3D = INTGRLCICO,X4D) 
X2D = INTGRL(ICO,X3D) 


XD = INTGRLC ICO, xX2D) 

OUTPUT = INTGRL(ICO,XD) 

FCNER = E**2 

El = INTGRLC ICO, FCNER) 
So ekOL FINTIM = 3.0 
TERMINAL 

COST = El 
END 
STOP 


FORTRAN 

* INITIALIZATION OF PARAMETERS FOR HOOKE SUBROUTINE 
TMPITECIT REAL*8(A-H,O-Z)} 
POPE NSION X(3);STEP(3) ,QC3),0003),WC3) 
eee!) = .001 

. 0001 

OOO | 


os. 


MINMAX = -1 

CALL HOOKE(X,STEP,N, ITMAX,CFTOL, ALPHA, BETA,CF ,Q,QQ,W, 
LIPRINT ,MINMAX) 

Ser 

END 
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Kivdkedvdkedetdeavaeavacdescacdedkeavdcakeakdesedeatdcskeskededcdescacakcakcacacakcacdkcdcacdvoedcdcakcdescacacsesescacsde sede ake 


em THESIS RESEARCH : 


FUNCTION MINIMIZATION 3RD ORD SYS 


* 


* THIS PROGRAM MINIMIZES THE COST FUNCTION, THE INTEGRAL* 


* OF ERROR SQUARED PLUS A WEIGHTING FACTOR (Q) TIMES 
* THE CONTROL EFFORT SQUARED AND RETURNS FEEDBACK 


* GAINS REQUIRED TO MINIMIZE THIS COST FUNTION. 


SiedeakeaedevesevededdesekdededeiskikkRseedelekickkkikdecdtedvkvleckkicctkkkkckeeKaRaK Te 


* 
D 


TITLE MINIMIZATION OF INTGRLCE**2+Q*U**2) ,3RD ORD SYSTEM 
oc 


COMMON/HANDJ/FLAG ,COST,KO,K1,K2 


mae) = KO, X(2) = K1, X(3) = K2 


INCON 


CONST K00=42.698049 , K10=5.2332549 , K20=. 30510201 


ICO=0. 


CONST Q = 0.00001 


INITIAL 
PFC FLAGAGE,.0.0)GO TO 5 
Kil = K10 
K2 = K20 
KO = KOO 
5 CONTINUE 
FLAG = FLAG + 1 
iy 
aku VATIVE 
R = 1.0 * STEP(0. 0) 
Poh OU LPT 
I = KO *E - K2 * X2D 
U= 1 - Ki*xp 
X3D = 96 * U -8*X2D - 12 *XD 
X2D = INTGRLCICO,xX3D) 


XD = INTGRL(CICO,X2D) 
OUTPUT = INTGRL(ICO,XD) 
FCNER = Eee 2+QeUr? 
E1 = INTGRLCICO,FCNER) 


CONTROL FINTIM = 3.00 
TERMINAL 


END 
STOP 


os 


COST = El 


FORTRAN 
* INITIALIZATION OF PARAMETERS FOR HOOKE SUBROUTINE 


PPLICIT REAL*8( A-H, 0-2) 
PENS VON <XC3),S5TEPC3) ,OC3) ,Q0(3) .WCs) 


STEP(1) = .001 
STEP(2) = .00001 
STEP(3) = 0.000001 
N = 3 

ITMAX = 500 

CFTOL = .00001 
APHA = 2.0 

Pena = 0.5 

IPRINT = 0 


se 


* 
hy 


MINMAX = -1 

CALL HOOKE(X,STEP,N, ITMAX,CFTOL, ALPHA, BETA,CF,Q,QQ,W, 
LIPRINT ,MINMAX) 

STOP 

END 


KKK kkk kee dekdkesdkedkdkdeakdkakdedcakakbdekakaesdkeskakedeakeskaedesakeadkeakese ke kdededededk eves sk 


* THESIS RESEARCH : FUNCTION MINIMIZATION 4TH ORD SYS * 
* THIS PROGRAM MINIMIZES THE COST FUNCTION, THE INTEGRAL* 
% OF ERROR SQUARED PLUS A WEIGHTING FACTOR TIMES THE % 
% CONTROL EFFORT SQUARED AND RETURNS THE FEEDBACK % 
* GAINS REQUIRED TO MINIMIZE THIS COST FUNTION. * 
Mretedevedetedede Re Le Dede se Ve Fe Se Ts MeN TEN TENS SENT TTF DETER TED TE DETER PORTER RETIRE R FRR 
* 

D COMMON/HANDJ/FLAG , COST ,KO,K1,K2,K3 

TITLE MINIMIZATION OF INTGRLCE**2+Q*U**2) ,4TH ORD SYSTEM 

a 


* X(1) = KO, X(€2) = K1, X(€3) = K2, X(€3) = K3 

CONST ICO=0. ,KO0=793. 83674 , K10= 87. 119566 ,K20=4. 0852083 
CONST K30 = .093384059 

CONST Q = 1. 0E-06 


INITIAL 
IF( FLAG. GE.0.0)GO TO 5 
KO = KOO 
Kl = K10 
K2 = K20 
K3 = K30 


5 CONTINUE 
FLAG = FLAG + 1 
% 


DERIVATIVE 


R = 1.0 * STEP(0.0) 

E = R - OUTPUT 

M1 = KO *E - K3 * X3D 

I = Ml - K2 * X2D 

(e= 1 = hieew XD 

X4D = U * 493.1 - 17 * X3D - 84 * X2D - 108 * XD 
X3D = INTGRL(ICO,X4D) 

pe INMIGR LC ICO oD) 


XD = INTGRLCICO,X2D) 
OUTPUT = INTGRLCICO,XD) 
FCNER = Es 2+Q"US% 2 
El = INTGRLCICO,FCNER) 
eenekOL FINTIM = 3.0 
TERMINAL 
COST = El 
END 
STOP 


ale 


FORTRAN 
* INITIALIZATION OF PARAMETERS FOR HOOKE SUBROUTINE 
PMPEICTT REAL*8( A-H,O-Z) 
DIMENSION X(4) ,STEP(4) ,Q(¢4) ,QQC4) ,W(4) 
O1 


Sebel) =". 
STEP(2) = .001 
SlEPCS) = 2001 
STEP(4) = .0001 
NS 


ITMAX = 500 


ee 


CFTOL = .000001 
ALPHA = 2.0 
BETA = 0.5 

BER INT “=20 
MINMAX = -1 


CALL HOOKE(X,STEP,N,ITMAX,CFTOL, ALPHA, BETA,CF ,Q,QQ,W, 
1IPRINT,MINMAX) 

STOP 

END 
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* THESTS RESEARCH: UNIT STEP RESPONSE 4TH ORDER SYSTEM % 
* THIS PROGRAM PRODUCES THE UNIT STEP RESPONSE FOR THE ¥ 
¥* MINIMIZED SYSTEM DESIGN AND THE BUTTERWORTH DESIGN. * 
Terese sede sess sevesede devs seve deve de Tete se wese se seve deveseve dese se ve sere Teese Nese He ee eck see sieve vest le sk teak seseskeaeae 
+) 
es 


eee SULT ERWORTH MINIMIZED SYSTEM UNIT STEP RESPONSE 
ae 


INCON ICO = 0.0 


PARAM KO = 99.482787 ,K1 =17.718436 ,K2 =1. 45109580 
PARAM K3 = .49327812E-01,WO = 14.89,Q = 0.0001 
DERIVATIVE 

R = 1. 0*STEP(0. 0) 

EB = R-BUT 

EM = R-MIN 


* 
* BUTTERWORTH DESIGN FROM STANDARD FORM 


X4B = WO**4*EB-2. 6*WO*X3B-3. 4*WO**2*X2B-2. OF WOK3*XB 
X3B = INTGRLC ICO, X4B) 

X2B = INTGRL(ICO,X3B) 

XB = INTGRLCICO,X2B) 

BUT = INTGRLC ICO, XB) 

EB1 = EB**2 


BCOST = INTGRL(ICO,EB1) 


ac 


*MINIMIZED DESIGN FROM FUNCTION MINIMIZATION 
Pile SEO = KS) ASD 
I = Ml - K2 * X2D 


U=T1 - K1*xXD 

X4D = 493.1 * U - 17 * X3D - 84 * X2D - 108 * XD 
X3D = INTGRLC ICO, X4D) 

X2D = INTGRLC ICO, X3D) 


XD = INTGRLC ICO, X2D) 
MIN = INTGRLCICO,xXD) 
EM] = EM**2 + Q*Us*2 
MCOST = INTGRL(CICO,EM1) 
CONTROL FINTIM = 3.00 
eee 01, TINE,BUT,MIN,BCOST,MCOST 
Poo. | ,BUT,MIN,MCOST,BCOST 
GRAPH (G1,DE=I1BM3279,LOGO=N) TIME(CUN=SEC) , BUT,MIN 
LABEL (G1)UNIT STEP RESPONSE 4TH ORDER SYSTEM ,Q =0.0001,WO=14. 89 
END 


PARAM WO=19. 85 ,KO0=309. 94665 ,K1=41. 110310, K2=2. 6116583 ,K3=. 072434062 
PARAM Q = 0.00001 

CANCEL LABEL(G1) 

LABEL(G1)UNIT STEP RESPONSE 4TH ORDER SYSTEM ,Q =0. 00001 ,WO=19. 85 
END 

rg 
PARAM WO=11. 15, K0=31. 353198 ,K1=7. 4873020 ,K2=. 77970130,K3=. 030946171 
PARAM Q = 0.001 

CANCEL LABEL(G1) 


a 


CABEL(G1)UNIT STEP RESPONSE 4TH ORDER SYSTEM ,Q =0.001,W0=11. 15 
END , 
STOP 


teksvavdeskoescavaeckscsvss deskatssckokseskseskakskskskcsvokoksksvckskskcsvckesakdcsdeskeskvoksvsesvslcakdeskslesksksesteskde ses se 


peelaeslS RESEARCH: UNIT STEP RESPONSE 4TH ORDER SYSTEM: * 
* THIS PROGRAM PRODUCES SATURATION RESPONSE FOR COMPARISON * 
* OF BODE AND MINIMIZED DESIGNS. * 


seveavacscskeakakseskvsksesvsestvakoktavsleskcskeskvsvseseskscvesteseskscscskclccstsvseskovsvsvseseseskcsesksksksletkvestesesksesestk se 


rig 
chy 
TITLE: SATURATION STUDY BODE,MINIMIZED SYSTEM DESIGNS,STEP RESPONSE 
r 


INCON ICO = 0.0 

ARRAY FN(4),FD(4),PN(1),PD(5) 

TABLE FN(1-4)=811. 64,13812. ,68248. ,87748 

TABLE FD(1-4)=1,240,8375, 75000 

TABLE PN(1)=493. 1 

TABLE PD(1-5)=1,17,84,108,0 

CONST WO=11. 15, KO=31. 353198,K1=7. 4873020 ,K2=. 77970130 ,K3=. 030946171 
CONST M=1.0 ,Q = .001 


’ 
riy 


DERIVATIVE 
R = 1. O*¥STEPCO. 0) 
EBA= R-BODEA 
* EM = R-MIN 
EMA = R-MINA 
¥ EMB = R-MINB 
* BODE DESIGN NO SATURATION 
# UB= TRNFR(3 ,3,0.0,FN,FD,EB) 
* BODE = TRNFR(O,4,0. 0,PN, PD, UB) 


* X1 = EB*2 + Q*UB2 

* Chee NGRLC TCO, x1) 

ee IEB = ERB**2 

ve IUB = UB***2 

2 IEEB INTGRICICO,1EB) 

* IUUB INTGRLC ICO, IUB) 

* BODE DESIGN WITH SATURATION 
UBA= TRNPR(3,3,0.0,FN,FD,EBA) 
LBUA = LIMIT(-2. 00,2. 00,M*UBA) 
BODEA = TRNFR(CO,4,0.0,PN,PD,LBUA) 


X1A = EBA‘*2 + Q&LBUAt*2 
CBA = INTGRL(ICO,X1A) 
IEB = EBA‘2 

IUB = LBUA***2 

IEEB = INTGRL(ICO,IEB) 
IUUB = INTGRL(ICO,IUB) 


ate 
#% 


Mee OD RESPONSE NO Se ern 
p= oti KO. =.K3 * X3D 

: i= ttl = K2.* XD 

- Ui = eS e Kk 1xD 


% K4D = 493.1 * UM- 17 * X3D - 84 * X2D ~- 108 * XD 
e X3D = INTGRLCICO,X4D) 

w X2D = INTGRLCICO,X3D) 

7 XD = INTGRLC ICO, X2D) 

¥ MIN = INTGRLCICO,XD) 

— E3 EM****2 +Q°*sUM***2 


Cie— aN TGRLCICO,ES) 
ReMi ess 2 


ri 1UM = UM*7%2 
IEEM = INTGRL(ICO,IEM) 
* IUUM = INTGRL(ICO, IUM) 


*MINIMIZED DESIGN WITH SATURATION 
M1A = EMA*KO - K3 * X3DA 
IA = MIA - K2 * X2DA 
UA = IA - K1*XDA 


LMUA = LIMIT(-2. 00,2. 00,M*UA) 

X4DA = 493.1 * LMUA - 17 * X3DA - 84 * X2DA - 108 * XDA 
X3DA = INTGRL(ICO,X4DA) 

X2DA = INTGRL(ICO,X3DA) 


XDA = INTGRL(ICO,X2DA) 
MINA = INTGRL(ICO,XDA) 


E3A = EMA**2 
CMA = INTGRL(ICO,E3A) 
IEM = EMA**2 
IUM = LMUA**2 
IEEM = INTGRL(ICO,IEM) 
IUUM = INTGRL(ICO,IUM) 


CONTROL FINTIM = 2.00 

SAVE .01,TIME,MINA,BODEA,CMA,CBA,LMUA,LBUA, IEEB , IUUB, IEEM, IUUM 
PRINT .1 ,MINA,BODEA,CMA,CBA, LMUA, LBUA, IEEB, IUUB, IEEM, IUUM 
GRAPH (G1,DE=1BM3279,LOGO=N) TIME(UN=SEC) ,MINA, BODEA 

LABEL (G1)STEP RESPONSE 4TH ORDER SYSTEM ,Q =0.001,WO=11. 15,L=2 


GRAPH (G2,DE=IBM3279,LOGO=N) TIMECUN=SEC) , LMUA,LBUA 

LABEL (G2)STEP RESPONSE 4TH ORDER SYSTEM ,Q =0.001,WO=11. 15,L=2 
END 

Sie 


deve duededcuedevssededevedededevede ve deve seve ve sek tested ve sede seve ve de Fe se Fe ede Fede Ke Fede Ke seve Fe seve sek Tek he seve seve 
* THESIS RESEARCH: UNIT STEP RESPONSE 3RD ORDER SYSTEM 7 
feeiiis PROGRAM PRODUCED THE UNIT STEP RESPONSE FOR ¥ 
* COMPARISON OF BUTTERWORTH AND MINIMIZED DESIGNS. ¥ 


sesevesededesescvededcdevevededesesesvevedsdesesesedesledevevesedededesedededesedevedevededededesevesksedevdesesicsesese sese ve 
* 


wc 
ieette: BULIERWORTH MINIMIZED SYSTEM STEP RESPONSE , COMPARISON 
* 


INCON ICO = 0.0 


PARAM WO=17. 12,K0=54. 906603 ,K1=5. 7501480 ,K2=. 23870600,Q = .0002 
PARAM WO=26. 72,K0=343. 75667 ,K1=18. 415644 ,K2=. 64320603,Q = .000005 
PARAM WO=15. 96,K0=42. 958667 ,K1=5. 2368362 ,K2=. 24983330,Q = .0005 
w 
DERIVATIVE 

R = 1. O¥STEPC(O. 0) 

PB = R=BUT 

Ebi Kooi 


*BUTTERWORTH DESIGN 
X3B = WO**3*EB-2. OFWO*XZB-2. OFWORH2*XB 
X2B = INTGRLC ICO,X3B) 
XB = INTGRLC ICO, X2B) 
BUT = INTGRLC ICO, XB) 
EB1 = EB*72 
BCOST = INTGRLCICO,EB1) 
ve 
*MINIMIZED DESIGN 
I = EM*KO - K2 * X2D 
U =I - K1*XD 
Pole 960 U- 8 * X2D - 12°" xD 
Fopr= INTCRIC ICO. X3D) 
XD = INTGRL(ICO,X2D) 
MIN = INTGRLCICO,XD) 
EM1 = EM*:52 + QU? 
MCOST = INTGRLCICO,EM1) 


’ 


iy 
CONTROL FINTIM = 2.00 

SAVE .01,TIME,BUT,MIN, BCOST, MCOST 

PRINT .1,BUT,MIN,BCOST,MCOST 

GRAPH (G1,DE=IBM3279,LOGO=N) TIME(UN=SEC) ,BUT,MIN 

LABEL(G1)UNIT STEP RESPONSE 3RD ORDER SYSTEM ,Q =0.0005,WO=15. 96 
END 

PARAM WO=17. 12,KO0=54. 906603,K1=5. 7501480,K2=. 23870600 

CANCEL LABEL(G1) 

LABEL(G1)UNIT STEP RESPONSE 3RD ORDER SYSTEM ,Q =0.0002,WO=17. 12 
END 


PARAM WO0=26. 72 ,K0=343. 75667 ,K1=18. 415644, K2=. 64320603 

CANCEL LABEL(G1) 

LABEL(G1)UNIT STEP RESPONSE 3RD ORDER SYSTEM ,Q =0. 000005 ,WO=26. 72 
END 

eOP 


$3 
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* THESIS RESEARCH: UNIT STEC RESPONSE S3RDMORD Gwe) > len ¥ 
%* THIS PROGRAM PROVIDES SATURATION RESPONSE FOR BODE AND ¥ 
* MINIMIZED DESIGNS FOR COMPARISON. * 
FHicdededevededetevevedededededevede seve destedede se vevevevewede ds vee seve seve te sed sesede seve de seston dese ve a seve seve seus 
¥ 
TITLE: BODE ,MINIMIZED SYSTEM SATURATION RESPONSE 3RD ORDER SYS. 

ye 

INCON ICO = 0.0 

CONST Q =0. 0005 

CONST KO = 42. 958667 ,K1=5. 2368362 ,K2 =. 24983330 ,WO=15.96,M=1.0 
ARRAY FN(3),FD(3),PN(1),PD(4) 

TABLE FNC Vee 6 — 104.17 3630-50. 1 2500e 

TABLE FD(1-3)=1.0,70.0,1000. 0 

TABLE PN(1) = 96 

TABLE PD(1-4)=1,8,12,0 


ac 


ae 


DERIVATIVE 
R = 1. O*STEPCO. 0) 
EB = R-BODE 
EM = R-MIN 

* EMA = R-MINA 

“ EMB = R-MINB 


*BODE DESIGN NO SATURATION 
UB= TRNFR(2,2,0.0,FN,FD,EB) 
BODE = TPRNERGO Ws 3020. N PDA WED 
X1 = EB? + Q*UB**2 
CB = INTGRLCICO,xX1) 

IEB = EB**2 

= [JB--<2 

TEES INTGRLOICOS TEED 
IUUB = INTGRL(ICO,IUB) 

* MINIMIZED DESIGN NO SATURATION 
PS Ee Oe ah ee 
UM= I - K1*XD 
X3D 96 * UM= 6G 32D ee  ae 
X2D INTCREGIGOR Mors 
XD/="INIGREC ICO ..on) 
MIN = INTGRLC ICO, XD) 

El = EM ***2+Q"UM*7*2 


CM INTGRL( ICO ,E1) 
IEM = EM**2 

IUM = UM***2 

IEEM INTGRL( ICO, IEM) 


IUUM = INTGRLC ICO, IUM) 


* BODE DESIGN WITH SATURATION 
+ UBA= TRNFR(2,2,0.0,FN,FD,EBA) 


ve LBUA.= LIMIT(=2700, 2.007 UBAY 

a BODEA = TRNFR(0,3,0.0,PN,PD,LBUA) 
% X1A = EBA**2 + Q*UBA**2 

* CBA = INTGRL(CICO,X1A) 


84 


* IEB = EBA**2 
oe IUB = LBUA***2 
ve TEEB = INIGRL( [C0. 1EB) 
we IUUB = INTGRLCICO,IUB) 


sk 


*MINIMIZED DESIGN WITH SATURATION 


ae TA=> EPMAKkKO = KZ * X2DA 

Ya UA = IA - Ki*XDA 

* LMUA = LIMIT(-2. 00,2. 00,M*UA) 
% X3DA = 96 * LMUA - 8 * X2DA - 12 * XDA 
v X2DA = INTGRLCICO,X3DA) 

Je XDA = INTGRLCICO,X2DA) 

oy MINA = INTGRLC ICO, XDA) 

ve E2 = EMA***2 

% CMA = INTGRUCICO,E2) 

“ IEM = EMA**2 

v IUM = LMUA**2 

ae IEEM = INTGRLCICO,IEM) 

sc [UUM = INTGRLC ICO, 1IUM) 


ac 

CONTROL FINTIM = 2.00 

SAVE .01,TIME,BODE, MIN,UB,UM,CB,CM,1IEEB,IUUB,IEEM,IUUM 

pieeeien, | ,BODE,MIN,UB,UM,CB,CM,IEEB,IUUB, JIEEM,1UUM 

GRAPH (G1,DE=1BM3279,LOGO=N) TIME(UN=SEC) ,BODE,MIN 

LABEL (G1)UNIT STEP RESPONSE 3RD ORDER SYSTEM ,Q =0. 0005 ,WO=15. 96, L=NO* 
Ey: 

GRAPH (G2,DE=IBM3279,LOGO=N) TIME(UN=SEC) ,UB,UM 

LABEL (G2)UNIT STEP RESPONSE 3RD ORDER SYSTEM ,Q =0. 0005 ,WO=15. 96,L=NO* 
END 

SOP 
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